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^ ■ Abstract 

. Gravitational waves are considered as metric perturbations about a curved background metric, 
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oo 
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rather than the flat Minkowski metric since several situations of physical interest can be discussed 
by this generalization. In this case, when the de Donder gauge is imposed, its preservation under 
infinitesimal spacetime diffeomorphisms is guaranteed if and only if the associated covector is ruled 
by a second-order hyperbolic operator which is the classical counterpart of the ghost operator in 
quantum gravity. In such a wave equation, the Ricci term has opposite sign with respect to the 
wave equation for Maxwell theory in the Lorenz gauge. We are, nevertheless, able to relate the 
solutions of the two problems, and the algorithm is applied to the case when the curved background 
geometry is the de Sitter spacetime. Such vector wave equations are studied in two different ways: 
i) an integral representation, ii) through a solution by factorization of the hyperbolic equation. The 
latter method is extended to the wave equation of metric perturbations in the de Sitter spacetime. 
This approach is a step towards a general discussion of gravitational waves in the de Sitter spacetime 
and might assume relevance in cosmology in order to study the stochastic background emerging 
from inflation. 
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I. INTRODUCTION 



One of the longstanding problems of modern gravitational physics is the detection of 
gravitational waves, for which the standard theoretical analysis relies upon the split of the 
space-time metric gab into "background plus perturbations" , i.e. 

gab = lab + hab, (l) 

where is the background Lorentzian metric, often taken to be of the Minkowski form 
rjab, while the symmetric tensor field hab describes perturbations about jab- However, the 
background jab needs not to be Minkowskian in several cases of physical interest, nor it 
has to be always a solution of the vacuum Einstein equations. As a consequence, we are 
therefore aiming to investigate in more detail what happens if the background space-time 
(M, Jab) has a non-vanishing Riemann curvature. 

This issue has to be seriously considered from an experimental point of view since the 
gravitational wave detectors of new generation are designed also to investigate strong-field 
regimes: this means that the physical situations, where only the standard Minkowski back- 
ground is taken into account, could be misleading in order to achieve self-consistent results. 

"erometers have been built in the United 



In particular, several ground-based laser inter 



States (LIGO) Europe (VIRGO and GEO) [2I, 3|, and Japan (TAMA) and are now 
in the data taking phase for frequency ranges about lO^^kHz. However, new advanced 
optical configurations allow to reach sensitivities slightly above and below the standard 
quantum limit for free test-particles, hence we are now approaching the epoch of second 



and third [6|] generation of gravitational wave detectors. This fact, in principle, allows to 
investigate wide ranges of frequencies where strong field regimes or alternative theories of 
gravity can be considered 3, Isl, • Besides, the laser interferometer space antenna (LISA) 
lO] (which is mainly devoted to work in the range 10"'^ ~ lO^^Hz) should fly within the 
next decade principally aimed at investigating the stochastic background of gravitational 
waves. At much lower frequencies {10^^"^ Hz), cosmic microwave background (CMB) probes, 
like the forthcoming PLANCK satellite, are designed to detect also gravitational waves by 
measuring the CMB polarization pA,] while millisecond pulsar timing can set interesting 
upper limits in the frequency range between 10~^ ~ 10~^ Hz [12]. At these frequencies, the 
large number of millisecond pulsars detectable by the square kilometer array would provide 
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a natural ensemble of clocks which can be used as multiple arms of a gravitational wave 
detector [l3 |. 

This forthcoming experimental situation is intriguing but deserves a serious theoretical 
analysis which cannot leave aside the rigorous investigation of strong-field regimes and the 
possibility that further polarization states of gravitational waves could come out in such 
regimes. For example, if one takes into account scalar-tensor theories of gravity [?] or 
higher-order theories 8|, scalar-massive gravitons should be considered. This implies that 
the standard approach where gravitational waves are assumed as small perturbations (com- 
ing only from Einstein's general relativity) on a Minkowski background could be totally 
insufficient. On the other hand, the existence of these further polarization modes could 
be a straightforward solution of the dark matter problem since massive gravitons could be 
testable cold dark matter candidates as discussed in 

In this paper, we want to face the issue of the rigorous formulation of gravitational wave 
problem in curved backgrounds. In particular, we want to perform a consistency analysis of 
gravitational waves in the de Sitter spacetime. Achieving solutions in this maximally sym- 
metric background could constitute the paradigm to investigate any curved spacetime by the 
same techniques and could have interesting cosmological applications if a conformal analysis 
is undertaken as, for example in [q], where it is shown how the amplitude of cosmological 
gravitational waves strictly depends on the cosmological background. 

It is straightforward to show that, in a covariant formulation, the supplementary condition 
for gravitational waves can be described by a functional acting on the space of symmetric 
rank-two tensors hab occurring in Eq. (1). For any choice of one gets a different 
realization of the invertible operator P^^'^'^ on metric perturbations. The basic equations of 
the theory read therefore as 

Pat^'hcd = 0, (2) 

^a{h) = 0, (3) 

where P^b'^'^ results from the expansion of the action functional to quadratic order in the 
metric perturbations. A deep link exists between classical and quantum theory, since in the 
latter, the one-loop analysis depends on the functional determinant of -P^f,'^"', after requiring 
that all metrics in Eq. (1) are positive-definite, i.e. Riemannian. Our analysis will instead 
be Lorentzian and classical. 
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The layout of the paper is the following. Section 2 studies the de Bonder choice for 
and its preservation, while Sees. 3 and 4 deal with massless Green functions in de Sitter 
spacetime. This is done because the problem of preserving Eq. (3) under infinitesimal 
diffeomorphisms leads precisely to vector wave equations. These are solved by an integral 
representation or by separation of variables. This analysis prepares the ground for studying 
the wave equation on metric perturbations itself through separation of variables, in Sec. 6. 
Concluding remarks and open problems are presented in Sec. 7, while relevant details are 
given in the Appendix. 



II. PRESERVATION OF THE DE BONDER SUPPLEMENTARY CONDITION 

Our first concern is how to implement in a consistent way the choice of supplementary 
condition. In general relativity, this is taken to be of the de Bonder type (below h = ■j'^'^hcd) 

Mh) = (^hab - l^abh^ , (4) 

if one wants to obtain the standard covariant wave operator on metric perturbations, where 

denotes covariant derivative with respect to the background metric 7^5. Under infinites- 
imal space-time diffeomorphisms, the metric perturbations suffer the variation (the round 
brackets denoting symmetrization) 

Shab = V(a (Pb), (5) 

where is a covector, with associated one-form ip^dx'' and vector field V'"^ (having set 
tpo. ^ 7"*(y9fe, which results from the isomorphism between tangent and cotangent space to 
the background space-time, that turns covectors into vectors, or the other way around). 
The change suffered from the de Bonder gauge in (4) when metric perturbations are varied 
according to (5) is then found to be 

5^a{h) ^ - {S^'U + Ra) fb, (6) 

where Q is the standard d'Alembert operator in curved space-time, i.e. Q = j'^'^'VcVd- 
By virtue of Eqs. (4) and (6), if the de Bonder gauge was originally satisfied, it is preserved 
under space-time diffeomorphisms if and only if (fh solves the equation 5^a{h) — 0. At this 
stage, to fully exploit what is known about the wave equation for Maxwell theory in curved 
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space-time in the Lorenz gauge^ we bear in mind that this reads as 



(7) 



This suggests adding iph to both sides of 5$a(^) = (see (6)), so as to cast it eventually 
in the form 

P," = 2R'^ y,,, (8) 

where 

Pa'^-^aU+Ra (9) 

is the standard gauge-field operator (see round brackets in Eq. (7)) in the Lorenz gauge. 
For this operator, the inverse is an integral operator with kernel given by the photon 
Green function, so that we can solve Eq. (8) in the form 



1°) + 2P " R, 



where (p^c^ is a solution of the homogeneous wave equation [17 1 



i'^Vf = 0. 



while is the inverse operator, satisfying 

pa pb ^^b 



(10) 



(11) 



(12) 



This is an integral operator with kernel given by the massless spin-1 Green function 
G'a6(.i;,x') = Gab'- The latter can be chosen, for example, to be of the Feynman type, i.e. 
that solution of the equation (see Appendix for the notation) 

5(x, x') 



(-^aD+^a) Gbe = ga 

having the asymptotic expansion as a ^ HQ 



Gab' 



gab' 



a + i£: 



+ Vab' log((T + is) + W, 



ab' 



where a{x, x') is the Ruse-Synge world function 
geodesic distance /i between the points x and x' . 



(13) 



(14) 



21 



22| , equal to half the square of the 



^ In [l6|. the author L. Lorenz, who was studying the identity of the vibrations of Ught with electrical 
currents, built a set of retarded potential for electrodynamics which, with hindsight, can be said to satisfy 
the gauge condition V^A^ = 0, which therefore should not be ascribed to H. Lorentz. 
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III. MASSLESS GREEN FUNCTIONS IN DE SITTER SPACETIME 



This general scheme can be completely implemented in the relevant case 23|] of de Sitter 
space where, relying upon the work in [2^1], we know that the massless spin-1 Green function 
reads as 

Gab' = <y{fi)gab' + (3{^)nany, (15) 

where x') = ^/2a{x, x') is the geodesic distance between x and x', n°(x, x') and n"-\x, x') 
are the unit tangents to the geodesic at x and x', respectively, for which 



na{x,x') = Va/i(x,x'), na'{x,x') = Va'/i(a;, x') , 



(16) 



while, in terms of the variable 



.^i(l + cos^,. 



the coefficient functions a and (3 are given, in four dimensions, by 24j ] 

1 



a{z) 



487r2p2 



+ - + - + ^ log(l - z) 



(3{z) 



1 



1 — z) ' z ' \z ' z'^ 

1 /I 1 



247r2p2 



1 



z \z z" 



log(l - z) 



(17) 



(18) 



(19) 



Strictly speaking, the formulae (18)-(19) are first derived in the Euclidean de Sitter space. 
In the Lorentzian de Sitter spacetime M which is what we are interested in, one can define 
the set [241 

Jj, = {x G M : 3 geodesic from x to x'} . (20) 

Moreover, it is well-known that M can be viewed as an hyperboloid imbedded in fiat space, 
i.e. as the set of points G R"^^ such that 



where rjab = diag(— 1, 1, 1), so that its induced metric reads as 



(21) 



ds' = r^abdVdY^ 



(22) 



As is stressed in Ref. 



24j |. the relation 



z{x^ x') 



r]abY''{x)Y^{x''' 



1 + 



(23) 
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is well defined both inside and outside Jx-, and it is an analytic function of the coordinates 
Y"'. Thus, Eq. (23) makes it possible to define z{x,x') everywhere on de Sitter, and one can 
define the geodesic distance 

/i(a;, x') = 2pcos^-^(V^) (24) 



as the limiting value 2S\ above the standard branch cut of cos~^. Along similar lines, the 
equations defining Wq, and Qav have right-hand sides which are analytic functions of the 
coordinates and are hence well defined everywhere on Lorentzian de Sitter spacetime 



24|. 



IV. EVALUATION OF THE KERNEL 



In a de Sitter background the Ricci tensor is proportional to the metric through the 
cosmological constant: Rab = ^gab, and hence the formulae (10), (15), (18) and (19) lead to 
the following explicit expression for the solution of the inhomogeneous wave equation (8): 

ip,{x) = ip^^\x) + 2A [ \a{z{fi{x, x')))gf + (3iz{fi{x, x')))n,n'''] ^a' {x') ^/-^ix')d^x' , (25) 



where, from Eq. (24), 



. n . -1 , VabY^ix)Y\x') 

/i(x,x ) = 2pcos ^2 



(26) 



while Eqs. (18) and (19) should be exploited to express a and P, bearing in mind Eq. (26) 
jointly with 



z{x,x) = - 



1 + COS 



p 



Moreover, the bivector g^' in the integrand (25) is given by 



24| 



(27) 



(28) 



p sin (/i/p)' 

The right-hand side of the formula expressing g^' is an analytic function of the coordinates 
Y"" and is therefore well defined everywhere on de Sitter 2J] . The integral on the right-hand 
side of Eq. (25) can be conveniently expressed the form 



fc{x) 



a(^)C-i(p)V,V"'p + (/5(^) -a(^))(V,p)(V'^'p)l^„.(a;')y^7M^^V, (29) 
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with a and /3 — a given by (cf. (18) and (19)) 



a{z) 



(1 + 2z) 
487rV 

1 



1 1 , . 

+ -log(l-2;) 



z{l-z) z^ 



3 + 2z-2z^) 3 

i^l°s(l--) 



487rV 

Equation (25) is therefore an integral equation reading as 
with unbounded kernel given by 



=2 



a(z)C-i(/i)VeV'^V + {P{z) - a(z))(V,/i)(V'^'/i) 



This kernel is indeed unbounded by virtue of the limits 



A^Ti'^ Wm. za{z) = -, 

2^0 / 



487rV lim(l - z)a{z) = 1, 

2^1 



487rV limz{(3{z) - a{z)) 

Jl „2 



3 
2' 



487r>^ lim(l - z){i5{z) - a{z)) = -3. 



At this stage, we can exploit (23) and (33) to re-express the kernel in the form 



, ^ (Vc^)(V°'^) 
247rV(l - z) 



z 



1 1 , , 



+ - log(l - z) 



(30) 
(31) 

(32) 



(33) 



(34) 
(35) 

(36) 
(37) 



(38) 



z{l — z) z^ 

Note now that ijf^ix) in Eq. (32), being a solution of the homogeneous vector wave equation 
(11), admits the Huygens' principle representation [181] 



(39) 



where, from Sec. 3 and the present section. 



Gcb';m' 2^-^cb';m'- 



(40) 
(41) 



Unlike the work in [18|, we here advocate the use of the Green function (15) rather than 
the sum, over all distinct geodesies between x and x' , of the Hadamard functions. To lowest 
order in the cosmological constant A, Eq. (39) may be used to approximate the desired 
solution of Eq. (32) in the form 

ip,{x) = <^(°)(a;) + A y K/VJ^ V^^tM^^V + OiK"). (42) 

Omitting indices for simplicity, the general algorithm for solving Eq. (32), here re- written 
in the form 

if = + A j Kip, (43) 

would be instead 

= <^(o)+A f Ky^m^ (44) 



y,„ = ^("^ + VA^ / ... / irV"\ (46) 



if = lim ipn. (47) 



V. SEPARATION OF VARIABLES FOR THE VECTOR WAVE EQUATIONS 

Consider now the de Sitter metric in standard spherical coordinates 

ds^ = -fdf + jdr^ + r\de^ + sin^ Odcj)^), (48) 

where / = 1 — iJ^r^ and H is the Hubble constant. This metric satisfies the vacuum Einstein 
equations with nonvanishing cosmological constant A such that = ^. Moreover, the 
timelike unit normal vector fields n to the t = constant hypersurfaces 

n^l (49) 

form a geodesic and irrotational congruence. The 3-metric induced on the t = constant 
hypersurfaces turns out to be conformally fiat, and the (1,1) form of the spacetime Ricci 
tensor is simply given by 

Ra = (50) 
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From the previous sections it is clear we are interested in the generahzed wave equation 



- UXa + ei?,% = + e 1^ Xa = 0, 



(51) 



where e = ±1. By virtue of the spherical symmetry of de Sitter spacetime thes e equat ions 
should be conveniently wtitten by using the expansion of X in vector harmonics 
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Following Zerilli 



we have 



X = Y{9)e 



-i{u)t—m4 



+ e 



+ ie 



—i{ujt—m4>) 



-i{Lut—m(p) 



'[foir)dt + Mr)dr] 

IDT 

dY 

-r sin6'/2(r)— + mf^{r)Y{9) 
dO 



(52) 



with Y{6) solution of the spherical harmonics equation 

.2 



m 



L 



Y = 0, 



(53) 



with L = + These equations lead to a system of coupled ordinary differential equations 
for the functions fo, fi, fs, besides a decoupled equation for /2 (/2 being related to the 
transverse part of X). Indeed, for I ^ 0,1 (the latter case being trivial), we have 



dr'^ 



r dr r"^ 



^V-L/ + 3/(l-/)(l-e) 



, 2z^(/-l) ^ 
Jo H 1. 111 



rf 



(54) 



d'h 2(2 - 3/) dfi 



dr'^ 



rf dr r'^f'i 



^V-L/-3e/(l -/) + /(/ -3) 



/i 



2MW) _ 2L 

XT iO r J 3i 

fj ^ J 



fh ^ 2(1 - 2/) df2 1 

dj.2 j.j dr 7-2y2 

d'fs 2(1 - /) df, 1 



c.V-L/ + /(l-/)(l-3e) /: 



^^r'_Lf + 3fil-f)il-e) /a - -/i- 

r 



(55) 
(56) 
(57) 



dr^ rf dr r'^f'^ 

Equation (56) for f^ir) can be easily integrated in terms of hypergeometric functions. In 



fact, assuming 

/2(r) = /-^-/(2^V(r), 
the resulting equation for ip reads as 



ar^ rj\ J dr f^j^ 



cuV^ - L + (1 - 3e)H\'^ + 3iujHr^ 



(58) 



(59) 
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the solution of which is, in general, of the form 



^(r) = Cir'sFi ( a_, a+, ^ + /, H\A + C^r-'-^F^ (a+, a_, ^ - /, i/V^ ) , (60) 



where we have defined 



— -3-2/±(13-12e)i/2 
H 



(61) 



Thus, when e = 1, which corresponds to studying the vector wave equation (7), one finds 

1 f2iuo 



a± = -- 



3 - 2/ ± 1 



4 V ^ 

whereas on taking e = — 1, i.e. our consistency Eq. (8), one gets 

1 f2iuj 



a-t = — - 



- 3 - 2/ ± 5 



(62) 



(63) 



4:\H 

The Lorenz gauge condition VqX" = 0, which only supplements Eq. (7), reduces instead to 



^ _r^fdh 2r(l - 2/) ^ 

J3 — T /i 



fo. 



(64) 



L dr L Lf 

The latter condition can be used, in principle, to obtain closed-form solutions of the various 
/o(r),/i(r),/3(r). 



VI. WAVE EQUATION FOR METRIC PERTURBATIONS 

Although the vector wave equations in de Sitter spacetime are already considerably in- 
volved, the final step consists in studying the invertible wave operator P^i^^'^ on metric per- 
turbations. On considering the DeWitt supermetric 



the de Bonder gauge in Eq. (4) can be re-expressed in the form 



(65) 



(66) 



and the resulting Lichnerowicz operator 29|], 30|| on metric perturbations, obtained by 
expansion of the Einstein-Hilbert action to quadratic order in hab, subject to $a(^) = 0, 



reads as 



3l| 



cl iz) d 



Pab"' ^ EJ\- □ + i?) - 2Ej^R\^f^"^ - Ej'R,^ - EJ^R, 
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(67) 



A wave equation for metric perturbations is therefore given by (see the Introduct 



where 



3l| 



D 2 



dn mc 
mil Inf 



Eab'" {-U + \R)-^Ejf5: 

Eab^' (- □ + \r^ + + |7a.7^' 

Eab^' {-U+\R)+^labr- 



The wave equation then becomes 



or 



implying also 



that is 



= Pab'^'hcd = + \r^ Kb + Jlabh, 

□ + -i? j /lab - -^labh = 0, 

-n\ + lRyh-lRh = o, 



after contraction with 7'^*. 
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A. Even metric perturbations 



Metric perturbations of even parity can be written in the form 



hi2 

^33 



fe-i(^t-r-'l^)Ho{r)Y{e), 



hoi 

ho3 = ime~'^''^-""f'^ho{r)Y{e), 
hii 



je-'^''*-'^'^^H2{r)Y{e), 

K{r)Y{e) + G{r 



dPY 



+G(r) 



-w?Y{e) + sin 61 cos 61 



dY 

lie 
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The wave equation (68) leads to the following system of coupled differential equations: 

l^ = 7f^'-'f '>'-^- ^^1'^''' + 2(1 - - 4/) - fL]H, 

2H\ dH2 dK 2H\L dG 
^ £ \ 

f dr f dr f dr 

2H^L AH^L 2H\l - 6/) ^ AH\ 

H 75 -"2 -^J^, 



2H^r dHo AH^r dK 2rLH^ dG 
-{ )L \ 

f dr f dr f dr 
| 2L(3-2/) ^ _ AL^^ _ 4(3-2/) ^^ _ 2H\l-2f) ^^^ 

^ = -^[.V - L/ - 4/(1 - /)].. - - Ih., 



1 + /,, 2 H^r 



H2 + —K + -—H^ 



rp rf p 



0, 



^-;^a-2/)f-;^l'^V-i^/ + 2/(3/-2)|G-^^ 

2 2H^ 
:H2 H ^Hf). 



r'f f 



To these equations one should add the de Bonder gauge components 



(75) 



2/ rV 2/ ^ / V 2 ; ' r/ 

^ 2(1-3/) ^ dHo ^ JK ^dG 2L^ 

dr rf dr dr dr r 

2L^ 4^^ 2H''r^^ 2iuj 

r^ r f f 

dh, 2(1-2/) L-2 ^^.^^irr rr. (^.. 

Metric perturbations of odd parity are instead found to vanish identically. 
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VII. CONCLUDING REMARKS 



A consistency analysis for gravitational waves in curved background is not, by itself, new 
in the literature and it has been outlined, for example, in Sec. II and Appendix B of Ref. 



32] • However, that paper was mainly concerned with gravitational instability in higher 



dimensions. In this paper, we have endeavoured to provide explicit solution formulae for 
the covector which solves the residual gauge problem expressed by Eq. (6), and this has 
been accomplished in Sees. 4 and 5. Section 6 has then moved on to work out all equations 
obeyed by metric perturbations subject to the de Bonder gauge in the de Sitter spacetime. 
Although the system (74)-(76) obeyed by metric perturbations looks very complicated, 
our approach is suitable for computer-assisted investigation of such equations, with the 
possibility of investigating gravitational waves in inflationary cosmology 



23| and/or ot 



ner 



curved backgrounds relevant for strong-gravity regimes. In particular, as discussed in [33| , a 



significant fraction of energy, in the form of a stochastic background of gravitational waves 
could emerge during the reheating after inflation. In this situation, the exact classification 
of gravitational wave solutions (as in Sec. 6) could be crucial in order to discriminate among 
the various signals. In a forthcoming paper, starting from the presented solutions, we will 
discuss the problem of generation and production of gravitational waves in the de Sitter 



background 



34|. 



APPENDIX A: BIVECTORS AND BISCALARS 



In Eq. (13), (7", is the geodesic parallel displacement bivector (in general, bitensors 
behave as a tensor both at x and at x') which effects parallel displacement of vectors along 
the geodesic from x' to x. In general, it is defined by the differential equations 



■,b „a ■,b' „a 

9c':b = ^ 9c':b' 



0, 



jointly with the coincidence limit 



lim gl, 



51. 



(Al) 



(A2) 



The bivector (7^,, when acting on a vector B'^' at x', gives therefore the vector which is 
obtained by parallel transport of 5^' to x along the geodesic connecting x and x' , i.e. 



9", i?" 



(A3) 
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In Eq. (14), A( biscalar built from the Van Vleck-Morette determinant 

D{x,x') = det{a,ab') (A4) 



according to 

A (x, x') = ^ D{x, x') (A5) 

The biscalar A{x,x') has unit coincidence limit: [A] = 1; as a function of x (resp. x'), it 

becomes infinite on any caustic formed by the geodesies emanating from x' (resp. x). When 



J]. 



A diverges in this way, x and x' are said to be conjugate points |3 
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